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SEARCH FOR CORRELATED UV AND X-RAY ABSORPTION 

IN NGC 3516 


Final Technical Report 
1. INTRODUCTION 

NGC 3516, a low-luminosity Seyfert galaxy, is one of a small fraction of Seyfert galaxies 
that exhibit broad absorption in a resonance line. Ulrich and Boisson (1983) first found 
that the C IV A1550, N V A1240, and Si IV A1440 emission lines all contain broad (EW 
~ 9A) absorption lines blueshifted with respect to the narrow lines. Voit, Shull, and 
Begelman (1987, hereafter VSB) used archival IUE data to demonstrate variability of the 
equivalent width, redshift, and line profiles on time scales as short as one week. In one 2 
week period, the line profile change had a P Cygni form, possibly indicating radial outflow 
of the UV absorbing gas. Since the absorption feature is often deeper than either the 
emission line or the continuum level alone, the absorbing gas must partially cover both the 
broad-line gas and the central continuum source. 

NGC 3516 also undergoes rapid X-ray spectral variations that have been attributed 
in part to variable photoelectric absorption (Halpern 1982; Reichert et al. 1985). Two 
models have been proposed to account for this behavior. In the “cold absorbers” model, 
clouds of moderate temperature and ionization parameter partially cover the X-ray source, 
with orbital motion producing the variable covering. The absorbers can be identified with 
the broad emission- line clouds, since the inferred column densities (~ 10 23 cm -2 ) and tem- 
peratures (T ~ 10 4 K) are comparable. An alternative model was proposed by Halpern 
(1984), invoking warm clouds, whose ionization balance can be strongly influenced by a 
variable nuclear continuum. When the ionization parameter is high, H and He are com- 
pletely stripped, and the soft X-rays are unattenuated, while there is significant absorption 
near 1 keV due to continuum absorption by 0 +6 and 0 +7 . Reichert, Mushotzky, and Holt 
(1986) proposed that the X-ray and UV absorption have a common origin in the broad 
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emission- line clouds. VSB argue against this hypothesis, based primarily on the large ve- 
locity coverage of the absorption line compared to that predicted in individual broad-line 
clouds. However, they do not rule out the possibility of common UV and X-ray absorption 
in gas other than the broad-line clouds. 

2. OBSERVATIONS 

In order to determine whether the UV and X-ray absorption in NGC 3156 are related, 
5 IUE observations were obtained in October, 1989, qusai-simultaneously with 4 Ginga 
observations. All the new IUE observations are listed in Table 1. The images were ob- 
tained with the large aperture in the low dispersion mode, resulting in a spectral resolution 
of ~ 5 A and ~ 8 A for the SWP and LWP, respectively. 

The spectra obtained on October 19 were severely contaminated by increased back- 
ground radiation due to solar flare activity. These spectra were found to be useless and 
were excluded from further data analysis. Four SWP and three LWP spectra of optimal 
quality were taken approximately equally spaced in time and spanning nine days all to- 
gether. In these spectra, we flagged the data points due to saturation, cosmic ray hits or 
the presence of reseaux. The coadded spectra are shown in Figure 1, and the the individual 
good SWP spectra are in Figure 2. Optimal exposure in the C IV A1550 line was achieved, 
and showed deep absorption. Only 3 out of the 4 Ginga observations were successful. 

These spectra were produced from the line-by-line images with the Gaussian extrac- 
tion (GEX) routines (Urry & Reichert 1988) at the Goddard Regional Data Analysis 
Facility. In short, these routines fit a Gaussian at each wavelength along the spatial di- 
rection of the image. Subsequently, these fits are required to vary smoothly along the 
spectral direction. This procedure minimizes the noise from the extraction process while 
maintaining an accurate flux level. At a distance of 27 Mpc (z = 0.009; H 0 = 100 km s -1 
Mpc -1 and ^ 0 =0 throughout this article), the spatial resolution of the SWP is 300 pc per 
line, which much larger than the bulk of the line emitting region in NGC 3516. Therefore, 
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the emission-line region, as well as the dominant nonstellar continuum, is detected as a 
point source and the use of a Gaussian fit is justified. 

Both the SWP and LWP data are corrected for the slight temperature dependency 
of the detector sensitivity. The LWP data were calibrated with the standards obtained in 
1984-1985 (Cassatella, Lloyd, & Riestra 1988). The May 1980 flux calibration was used 
for the SWP data. The accuracy of this calibration is estimated at ~10% (Bohlin & Holm 
1980), but work by Finley et al. (1990) suggests that the actual calibration could be up 
to 20% different from the May 1980 calibration. 

The photoelectric data in the Fine Error Sensor (FES) on board the IUE were con- 
verted to V-band magnitudes and listed in Table 1. The formal statistical error due to 
the calibration is ±0.07 mag (Sonnebom et al. 1987). The field of view of the FES has a 
diameter of \2" while NGC 3516 extends 2/1 by 1/8, therefore only the central emission 
region is detected. 


3. RESULTS 

3.1. The Short-term UV Variability 

From these spectra, we calculated the average fluxes in the continuum intervals and 
estimated the errors by using the standard deviation of the data points in the continuum 
interval. These errors might underestimate the true error, as the uncertainty in the cal- 
ibration is not taken into account. The continuum was studied for variability by fitting 
the data points with the average flux (shown in Figure 3 as the dashed lines). A change 
in continuum level is obvious between October 7 and 9 (day 281 and 283, respectively). 
Averaged over wavelength the change is 16% (the ratio varies from 88% ± 8% around 1710 
A to 81% ± 5% around 1955 A). The statistics of the test for variability on the whole data 
set are shown in Table 2, which indicates significant variability at the 99.6% level in the 
most sensitive (1940-1970 A) part of the SWP. 
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The variability comes from the Oct. 7 data, since exclusion of this day results in no 
significant variability (Table 3), with goodness-of fit parameter q close to 1. An upper limit 
on the variability excluding Oct. 7, conservatively estimated from the lcr error estimates 
on the fluxes, is found to be ±6% around 1955 A and ±7% around 3000 A on a timescale 
from days up to a week. The power law index of the combined SWP and LWP continuum 
showed no significant variations from the mean either. 

On average, the V magnitude was 12.31 on October 7 and 12.41 between October 9 
and 16; a small initial decrease and a constant optical brightness afterwards, consistent 
with the decrease detected in the UV continuum between October 7 and 9 and the steady 
ultraviolet flux level between October 9 and 16. The decrease in the UV is about 0.2 mag, 
while the optical decrease is only 0.1 mag. 

3.2. Emission Lines 

We coadded the spectra to study the emission lines. The measured line properties are 
listed in Table 4 and the lines are identified in Figures 1. The Lya emission line could not 
be measured, as the blue wing is blended with the geocoronal Lyct line and the red wing 
with the strong N V absorption and possibly the N V emission as well. Following Ulrich 
and Boisson (1983) we measured the intensities of the emission features, uncorrected for 
interstellar absorption, absorption lines (also for C IV) or line blends. In Table 4 we list 
the formal errors, where the error (relative to the flux) in the individual line emission data 
point is approximated by the relative <7 in the continuum bin bordering the emission line. 
The major source of error in the line measurements is the uncertainty in the continuum 
level, which is contaminated by the broad wings of the emission lines, Galactic absorption 
lines, and blended Fe II line and Balmer continuum emission. Therefore, the formal error 
seriously underestimates the actual error in the line measurements. 

These 1989 UV line measurements can be compared to the ones reported by Ulrich 
and Boisson (1983) obtained from 1981 spectra. Very good agreement was found for the 


-4 - 


strong, low excitation lines (C III], Si III], and Mg II) and the weak lines (N IV, 0 III], and 
N III]). However, the reported 0 III A3133 line was not found; we established an upper 
limit at 25% of the flux listed by Ulrich and Boisson. Also, different results were found 
for the strong, high excitation lines (Si IV, 0 IV, and C IV): these line fluxes were much 
lower, which could be accounted for by the significantly stronger absorption features found 
in these emission lines. Compared to the 1981 spectrum the absorption is much stronger 
in the C IV, N V, and Si IV, 0 IV lines. 

The ratio of C III]/C IV measured in NGC 3516 is 0.19. This is a very conservative 
upper limit as there is flux from Si III] included in the C III] measurement and the C IV 
is severely underestimated due to the strong absorption. Such a low ratio is exceptional 
for quasars or Seyfert 1 galaxies (Kwan and Krolik 1981) and as this ratio is used to set 
an upper limit on the density in the BLR. A ratio <C 0.19 would allow for high densities 
(^e ~ 10 10 cm 3 ). The strength of Si III] relative to C III] is expected to increase with 
higher densities as well (Rees, Netzer, & Ferland 1989). 

3.3. Galactic Absorption Lines 

We measured the absorption lines listed in Table 5 and indicated by tick marks in 
Figure 1. All these lines arise from low-ionization species. There are also strong absorption 
features close to the redshifted emission lines of N V and C IV. These absorption lines are 
associated with material within the NGC 3516 system. Considering the wavelength of the 
Si IV AA1393,1402 absorption doublet (A a b s are 1403il and 1412i2 A) This feature is 
also thought to arise in the Seyfert and not in our Galaxy. In contrast, the Mg II line is 
centered very close to the rest wavelength and therefore attributed to Galactic absorption. 

3.4. The C IV, N V, and Si IV Absorption Features 

We measured the equivalent width of the CIV, N V, and Si IV absorption lines and 
listed the results in Table 6. The errors listed for the line measurements are again formal, 
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and therefore lower limits to the true errors. Following VSB, the equivalent width of CIV 
was measured in two ways: first by joining the two peaks by a straight line and second, 
by fitting a Gaussian profile to the narrow emission line components (as the narrow-line 
profile in QSOs is approximated well by a Gaussian). The first method undoubtedly 
underestimates the true equivalent width, while the second method is highly uncertain 
in its speculative fit of the emission profile. The Gaussian fit peaks around 1555 A and 
not around the expected emission peak of 1564 A. There is evidence that the base of 
the emission feature is also blueshifted with respect to the systemic velocity and that the 
absorption feature is not shifted considerably compared to the C IV emission line profile 
(VSB). The N V line was fit by a horizontal line to its red emission peak and the Si IV 
by a straight line to its red and blue peak. The Si IV line is weak and is hard to measure 
accurately in the relatively strong noise. 

The overall C IV line profile is remarkably symmetric, which is exceptional as most 
C IV profiles in NGC 3516 are distinctly asymmetric with a stronger red than blue peak 
(VSB). Compared to the observations from 1978 to 1984, the 1989 continuum is at a 
high level, while the peaks of the C IV line reach the highest flux levels ever. The C 
IV absorption equivalent width is also very strong, measured either way, and the N V 
equivalent width is unsurpassed. 

There is an absence of significant variability in the October 1989 absorption measure- 
ments (Table 6). The height of the peaks deviated only for the blue peak on October 
9. This also accounts for the smaller equivalent width under the straight line fit found 
on that date. The equivalent widths of C IV (under the Gaussian fit) and N V showed 
no variation within their formal errors, thus certainly excluding significant variability. In 
this light, we are justified to coadd the spectra and study the absorption features in the 
coadded spectrum. The coadded spectrum has a total exposure time of 1260 minutes (far 
more than any earlier coadded spectrum studied; VSB, Walter et al. 1990), resulting in 
spectrum with an unequalled signal-to-noise ratio. 
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3.5. Lower Limit on the Carbon Column Density 

We measured the optical depth r derived from the ratio of the observed flux to the 
estimated, unabsorbed flux: 

r = -ln(F*/F“ nabs ). 

Again we approximate the unabsorbed profile with a straight line and a Gaussian fit to 
the narrow line component. This calculation is only correct (Barlow, Junkkarinen, and 
Burbidge 1989) if there is no scattered light and emission from the absorption line region; 
if the absorption line does not contain any unresolved narrow fines; and if the absorber 
covers the emitting source. The spectral resolution of the SWP camera is low, so any 
unresolved narrow fines cannot be excluded with certainty. VSB showed that the absorbing 
region must cover (at least partially) the fine emitting region and the continuum. Together 
with the uncertainty of the unabsorbed fine profile, this estimate of the optical depth is 
therefore slightly uncertain. The optical depth reaches values slightly greater than 1 in 
the line center. C IV absorption ranges from ~ 0 km s -1 to ~ —3000 km s -1 . The 
absorption peaks around —1000 km s L The C IV absorption falls off more rapidly on 
the red than the blue side and there is evidence (e.g., on October 16) that the absorption 
profile flattens out around the minimum. The N V absorption also sets in around the 
redshifted rest wavelength and reaches a minimum around -1000 km s -1 , but extends to 
at most ~ —3000 km s _1 . 

The column density of C 3 "*", Nq3+ , can be estimated from the optical depth (Lucy 
1971) using the method by Junkkarinen, Burbidge, and Smith (1983): 

j Td ^ x )- 

where u 0 is the rest frequency, and / the oscillator strength. Considering the uncertainty 
in r we can only determine a lower limit to N c s+. It is at least 2.2 x 10 15 cm -2 for the 
Gaussian fit, and 1.8 x 10 15 cm 2 for the straight fine fit. Such Nq 3 + would correspond 
to iV H > 2 x 10 19 cm -2 , assuming cosmic abundances and that 1/3 of all carbon is three 
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times ionized. VSB estimated a similar lower limit to N c s+ of ~ 10 15 cm -2 using the 
curve of growth method. 

3.6. Estimate of the Distance from the Absorber to the Continuum Source 

Continuing on the assumption that photoionization is the main ionization mechanism 
in the absorbing region, we can relate the ionization parameter U to n e , the density of the 
absorber, r, the distance from the absorber to the continuum source, in addition to observ- 
able quantities. For U we use the value at which the fractional abundance of C 3+ peaks, 
Up, as large deviations from this value are not expected considering the obvious strength 
of the C IV features. Smith and Penston (1988) calculated U p with the photoionization 
code CLOUDY and a = 1.5. For log TVh = 20, 21, and 22, they find log Up = -2.18, —1.65, 
and -0.83, respectively. This would give us in the case of NGC 3516 and log Nn = 20, 
n e r 2 ~ 2 x 10 43 cm -1 , or with n in 10 5 cm -3 and r in pc: r 2 c ~ 20 /n 5 . For log iV H = 21 
and 22, r 2 c ~ 6/n 5 and l/n 5 , respectively. 

3.7. Variability in the Continuum and Absorption 

NGC 3516 is variable, both in the continuum and absorption. Such variability has 
been reported for other Seyferts as well; in particular, NGC 4151 shows absorption in 
C IV which changes in equivalent width and wavelength with time (Clavel et al. 1987). 
Variability can be attributed to either ionization changes in the absorber or changes in 
the bulk motion of the absorber in the line of sight. Boroson et al. (1991) reported on 
the appearance of a new absorption system in Mkn 231, which these authors explained by 
motion of an absorbing cloud into the line of sight. Mkn 231 differs from NGC 3516 in the 
fact that it is not as highly ionized. We will show that the variability of NGC 3516 can 
readily be explained to be due to changes in the ionization of the absorber. 

We retrieved all the IUE archival data (up to 1989) and extracted these with the 
Gaussian routines as described above. The data obtained with the old software (before 
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1980) were reprocessed with the corrected ITF. We searched the archival data for vari- 
ability on the shortest timescale. To follow the actual change in the continuum and C IV 
absorption, observations taken frequently are required. Observations taken in 1984 showed 
the fastest change in the shortest period (see also VSB and Walter et al. 1990). Definite 
variability in both the continuum flux level and equivalent width of the C IV absorption 
was found on a time scale of about a week. In this episode, anticorrelation of the contin- 
uum flux and C IV equivalent width, as reported by VSB, was seen. Between October and 
November, the continuum flux decreased, while the absorption in C IV became stronger. 

In theory, if the absorbing region is in photoionization equilibrium, then the fractional 
abundance of C 3+ first increases with U, reaches a peak value, U p , and then decreases with 
the ionization parameter. In 1984, the continuum flux decreased, resulting in a smaller 
U. If we assume that the changes in the absorbing strength result from changes in the 
fractional abundance (which, in turn, is caused by the time variability of the continuum 
source), then an increase in the C IV equivalent width would correspond to an increase 
in the fractional abundance of C 3+ . This, together with the decreasing U implies that 
U > U p . The lower limit to U corresponds to an upper limit on n 5 Tp C . VSB inferred a 
lower limit for the density from the observed upper limit on the recombination time t rec 
during the 1984 epoch: n 5 ~ 10/f rec ~ 1 (with t iec in days). Together with r pc < (20 /n 5 ) 
derived above, this would lead to an upper limit to r of 5 pc assuming log N#= 20 (for 
log JVh = 21 and 22, r would be < 3 and < 1 pc, respectively). 

Smith & Penston (1988) also predicted the correlation properties between the absorp- 
tion from C IV, N V, and Si IV. Roughly, the absorption strength of N V and Si IV vary 
with C IV (VSB). However, it is difficult to assess whether there is a significant correlation 
with the C IV absorption: the Si IV feature is very weak, but appears to be correlated, 
while the stronger N V line deviates from a possible correlation during several epochs (see 
Figure 9 of Walter et al.). This is unfortunate, as the correlations and anti-correlations 
could constrain the ionization parameter further. 
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4. DISCUSSION 


4.1. A Comparison with BAL QSOs 

There is a class of QSOs, the Broad Absorption Line QSOs, which has much in 
common with NGC 3516 and could very well be considered to be its high luminosity, high 
redshift equivalent. Since most observed BAL QSOs have large redshifts (z > 1), they have 
been studied extensively at high resolution (~ 1 A) with optical telescopes. The BAL 
profile can be highly structured with multiple troughs. In comparison, the nearby galaxies 
exhibiting such absorbing features have been observed with IUE and display smooth 
absorption profiles, but these could be attributed to the low resolution and a structured 
profile cannot be excluded. 

The most obvious similarity between NGC 3516 and the BAL QSOs is the absorbing 
feature in the high ionization lines: C IV, N V, and Si IV. There are BAL QSOs with 
absorption in low-ionization lines (e.g., Mg II), but they are the exception to the rule. The 
most striking difference is the strength of the absorption line: the equivalent width is ~40 
A in the BAL QSOs and at most 10 A in NGC 3516. The range in velocity is also quite 
different: from ~ 0 to at least ~ 10, 000 km s -1 in BAL QSOs and from ~ 0 to ~ 3500 
km s -1 in NGC 3516. Variability has been reported for both NGC 3516 and several BAL 
QSOs, only the timescale differs considerably: NGC 3516 varies on a timescale of months 
if not weeks, while BAL QSOs have been reported to vary on a timescale of years. 

Based on the failure to detect broad [O III] emission (which is expected to arise in 
the BAL region through collisional excitation) in the nearby BAL QSO PG 1700+518, 
Tumshek et al. (1985) constrain n e to be larger than the critical density for de-excitation 
(6 x 10 5 cm -3 ), unless the covering fraction of the absorber is very small (q <0.05, which 
is not in agreement with the assumption that all QSOs have BAL regions). With a similar 
argument as presented above, this implies that the absorbing gas is located around r < 30 
pc in PG 1700. Assuming that the physical conditions are similar in NGC 3516 and 
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making use of the fact that this limit scales with the square root of the luminosity, this 
would correspond to an upper limit for NGC 3516 comparable to the one derived above: 

titoo < (6.3 x 10 -3 ) 1/2 x 30 < 2 pc, 

where the subscripts 3516 and 1700 refer to NGC 3516 and PG 1700+518, respectively). 

A lower limit to r in BAL QSOs is derived from the fact that Lya is occulted by the N 
V absorption, implying that the BAL clouds are outside the Lyo emitting region, i.e. r > 1 
pc. A similar argument can be made for NGC 3516, where the absorbers have to be located 
outside the region where the C IV and broad Ly a emission originates (but not necessarily 
outside the C III] and Mg II emitting region). In the case of NGC 4151 there is evidence 
for a stratified emission region, with the core of the C III] emission line originating far 
outside the C IV line region (Clavel et al. 1990). In the case of NGC 3516, similar studies 
of the emission line variations (Edelson, Krolik, and Pike 1990) did not give significant 
estimates on the emission line region, let alone on its stratification. However, this could 
be attributed to the uneven and undersampling of the observations of NGC 3516, and it is 
not unrealistic to assume a stratified emission region for NGC 3516. This would put the 
absorber outside the C IV emitting region, i.e., r >25 It-days (0.02 pc). NGC 4151 has a 
similar luminosity as NGC 3516 and therefore the estimates for the emission region derived 
in NGC 4151 can be applied directly to NGC 3516. Recent high-rate reverberation studies 
of NGC 4151, both optically and with IUE, also showed that the BLR is much smaller 
than previously thought: r < 21 lt-days for the Balmer lines (Maoz et al. 1991) and r < 7 
lt-days for the C IV line (Clavel et al. 1990). This will decrease the lower limit on r even 
further to > 0.005 pc. 

Next we address the fact that the absorbing features are weaker and narrower in 
NGC 3516 compared to the BAL QSOs. In a first approximation, this can be attributed 
to the much lower luminosity in NGC 3516 which cannot accelerate the absorber to the 
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high velocities seen in BAL QSOs, resulting therefore in absorption features with smaller 
equivalent widths and narrower breadths. We can envision two simple models: 

1. The radiative acceleration is the dominant force acting on the absorbing clouds. The 
terminal velocity v of the absorbing clouds which are accelerated from a starting distance 
r min scales as (L/r m in) 1 / 2 (Weymann, Turnshek, &; Christiansen 1985, hereafter WTC). 
For BAL QSOs J'min ~ 10 pc for intermediate velocities and r m j n ~ 1 pc for the high- 
velocity clouds. In particular, in the case of QSO 1303+308 (Foltz et al. 1987) a similar 
argument led to an estimate of r m j n ~ 3 pc for the high velocity clouds (v ~ 10,000 km 
s -1 ). This would correspond to a terminal velocity in NGC 3516 of 

1/2 

UQ 1303 ~ (3200) -1 / 2 x (3/0.01) 1 / 2 x 10, 000 ~ 3000 kins’ 1 . 

Associating the clouds, which have been accelerated from the minimum distance and have 
reached the maximum velocity, with the blue edge of the absorbing feature, we find terminal 
velocities of at least 3000 km s -1 in NGC 3516, in good agreement with the scaled down 
version of a BAL QSO. 

2. The acceleration is due to a wind, in which ram pressure is the dominant force. WTC 
noted that the observed profile of the BAL, with a sharp red edge and a slowly rising 
blue edge, could easily be interpreted as arising from the line of sight passing through the 
main clouds in addition to smaller clouds which have been torn off the main clouds and 
accelerated to higher velocities. NGC 3516 shows a similar profile in the sense that the red 
edge is steeper than the blue edge (Figures 1 and 2). The terminal velocity scales similar 
to the radiative acceleration, but the velocities are expected (for similar L and r m j n ) to 
be larger by a factor ~10 (WTC). Thus, the clouds are more efficiently accelerated by 
a wind than by radiation and therefore the clouds in the wind model can originate at 
larger distances. For this reason, and for the need to invoke a confining medium for the 
radiatively accelerated clouds, the wind model is favoured in BAL QSOs. The wind model 
is equally applicable in NGC 3516, allowing for the high velocity clouds to originate at 
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moderate distances and the intermediate velocity clouds (the majority as the absorption 
clearly peaks around ~1000 km s -1 ) to originate closer to the continuum source. The same 
problems which plague these acceleration models in BAL QSOs will have to be addressed 
for NGC 3516, e.g., the existence of the wind, and the instabilities of the clouds (WTC). 
VSB found evidence for the presence of a wind in NGC 3516: the C IV profile change in a 
two week period displayed a P Cygni form, indicating radial outflow of the UV absorbing 
gas. Also, Goad & Gallagher (1987) found optical evidence for outflow with velocities of 
~ 700 km s -1 , on a scale of 100 pc. 

4.2. The X-ray - UV Connection 

Since the X-ray continuum source is smaller than the UV source, and as the UV 
absorbers are (at least partially) covering the continuum, the material producing the UV 
absorption could be expected to produce the X-ray absorption as well. If the UV and X-ray 
absorption had been observed to be well correlated, then there would have been strong 
evidence that one component is responsible for both types of absorption. Unfortunately, 
neither the UV nor the X-ray absorption was found to vary during our simultaneous 
observations in 1989 and this correlation could not be established. 

There are at least a dozen AGN which have been observed with IUE and which display 
UV absorption features of intrinsic origin (i.e., not of intergalactic or Galactic origin; Ulrich 
1988, Kinney et al. 1990). Some are genuine BAL QSOs, others are nearby galaxies very 
similar to NGC 3516. Most have strong and variable X-ray emission ( A r j j ~ 10 22 cm -2 ), 
but there are galaxies with no larger Nr than expected for Galactic absorption (e.g., 
MCG 8-11-11). If the column density is close to the lower limit derived for NGC 3516, 
i.e. Nh ~ 10 20 cm -2 , observable UV absorption lines could be produced, but the X-ray 
absorption would be difficult to detect. Therefore, the presence of UV absorption lines 
together with the absence of substantial X-ray absorption, does not rule out that the UV 
absorption arises in the same location as the X-ray absorption (e.g., in the broad line region 
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as suggested by Reichert, Mushotzky, and Holt 1986). Nevertheless, there are highly X-ray 
absorbed galaxies which do not show any evidence for absorption in the UV. A complete 
study of the X-ray and UV absorption properties galaxies like NGC 3516 and BAL QSOs 
has not yet been undertaken. Therefore, it is too early to conclude from the absorption 
properties of AGNs in the UV and X-ray that these arise in the same location and through 
the same process. 

In our Ginga spectra of NGC 3516, as in other observations of Seyfert 1 galaxies 
(Piro, Yamauchi, & Matsuoka 1990), we found evidence for two absorption components: 
the ‘standard’ cool absorption with Ah ~ 10 22 ' 6 cm -2 , and the very thick cool absorber 
with Nft ~ 10 24 ' 6 cm -2 . Alternatively, the second component can be interpreted to arise 
through reflection from very thick cool matter (Ah > 10 25 cm -2 ; Lightman & White 1988). 
Either way it suggests the presence of dense material in the Seyfert galaxy. If the UV and 
X-ray absorption have the same origin, with which X-ray component is the UV absorption 
associated? Obviously, if the hard tail above 7 keV arises in a reflective slab, no absorption 
is expected, but the X-ray observations do not allow us to distinguish between this model 
and the very thick cool absorber alternative. Assuming that two X-ray absorbers are 
present, a thin and thick one, which could be responsible for the UV absorption? The IUE 
observations established a lower limit on Ah of ~ 10 20 cm -2 , allowing both components. 

We will show that the thick X-ray absorption could arise in a small number (~ 30) 
of dense clouds, which would also be responsible for the UV absorption. Invoked by the 
small BLR sizes found recently through reverberation studies, high-density models for 
the broad line clouds have been calculated (Rees, Netzer, & Ferland 1989; Ferland and 
Persson 1989), which show that for densities around 10 11 — 10 12 cm -3 , C IV would be 
stronger than thought previously for the standard, low density BLR models. The same is 
true for Si IV and NV, the other lines which are seen in absorption. The lack of Mg II 
A 2798 absorption is more difficult to explain. Its strength peaks around 10 10 cm -3 and 
levels off at higher densities faster than C IV, but still remains at considerable intensity. 
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Specifically, comparing the relative line intensities at 10 10 cm -3 and 10 12 cm -3 (for model 
A of Rees, Netzer, & Ferland 1989), the C III] line decreases by 90%, C IV increases 
by 7%, N V increases by a factor of 7, and Mg II decreases by 45%. If the absorption 
regions are highly ionized, little Mg+ would be expected, while C 3+ and N 4+ would be 
abundant. This is illustrated by the model in which the ionization parameter is allowed to 
vary (model D, Figure 3 of Rees, Netzer, & Ferland 1989). Mg+ and C 2+ diminish quickly 
with increasing ionization, while Si 3+ and C 3+ decrease slowly. Thus, the absence of the 
low-ionization species can be used to constrain the ionization parameter U. In the model 
of Smith & Penston (1988), this implies U > 0.006 and U > 0.05 for log Nn = 21 and 22, 
respectively. Our estimates of the lower limit to U fulfill this inequality (§3.7). Some care 
has to be taken with the interpretation of the emission line models to predict the absorption 
line properties. While there can be little doubt that the high-ionization species (like C 3+ ) 
are present in such high-density broad line clouds, these might be expected at the side 
illuminated by the continuum source. In order for absorption to take place in the same 
broad line clouds, the absorber has to be located further away from the high-ionization 
front. It might be difficult to achieve the right ionization level there as the photoionization 
flux will be absorbed by the intervening dense material. The same argument holds if the 
UV absorption is only due to the outer broad line clouds, as the intervening clouds will 
intercept the ionizing flux. 

As pointed out by VSB, in the standard broad-line cloud model at least 30 clouds 
are needed to achieve the observed width in the absorption lines, as the thermal velocity 
dispersion in a single cloud is small. The Nu of a standard dense cloud is ~ 10 23 cm -2 , 
therefore 30 clouds will result in the observed X-ray column density of 10 24-6 cm -2 . These 
clouds will have to be located quite close to continuum source. The ionization parameter 
increases with JVh and if we (as a rough estimate) extrapolate this trend as calculated 
by Smith V Penston (1988) to log A r n = 24, then log U p ~ 0.4. Following the same 
argument as above, this gives r <0.3 pc. Walter et al. (1990) argue that the UV absorber 


- 15 - 



is located in a single cloud. This will be difficult to achieve: in order to obtain the required 
width in the absorption line from one cloud, the internal velocity dispersion has to be 
enormous. Such a cloud would not be stable. Walter et al. arrived at the single cloud 
model from the fact that nine consecutive observations followed two (anti)correlations 
between the absorption strength and the continuum flux. However, seven observations 
do not follow these correlations. It is therefore difficult to maintain that the trend of 
decreasing absorption strength with increasing continuum flux has to be explained with 
the presence of a single cloud. 


5. CONCLUSIONS 

In conclusion, NGC 3516 has the characteristics of a low luminosity BAL QSO. The 
density of the absorbing region has a similar lower limit in this Seyfert as in the BAL 
QSOs. The absorber of C IV, N V, and Si IV is probably located between 0.01 and 5 pc 
in NGC 3516, much closer to the nucleus than in the high luminosity BAL QSOs. The 
clouds in NGC 3516 do not reach the large velocities and distances from the continuum 
source as in BAL QSOs, because the central engine in NGC 3516 can only accelerate them 
to moderate velocities. Therefore, the range of velocities and the maximum velocities are 
smaller, resulting in an absorbing feature which is narrower with a smaller blueshift in the 
rest frame of the Seyfert. 

Our simultaneous observations in the UV and X-ray are not in disagreement with 
the assumption that the UV absorber is also responsible for the X-ray absorption. This 
absorption could be identified with the dense component (ATh ~ 10 24 ' 6 cm -2 ), detected 
in the Ginga observations. If so, then the absorber would be located very close to the 
continuum source: r < 0.3 pc. The lack of Mg II absorption could be attributed to the high 
ionization level of the absorbing region. 
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6. CONTINUING AND FUTURE WORK 


There are two uncertainties which limit the usefulness of the current measurements of 
absorbing columns of C IV, N V, and Si IV. One is the shape of the unabsorbed continuum 
plus line emission. The other is the profile of the absorption features. The latter can be 
studied better with high resolution observations. Such observations for high 2 BAL QSOs 
showed a structured, detailed absorption through, in contrast to the profiles obtained with 
IUE. Such BAL profiles have been interpreted to arise from a large (possibly up to ~ 10 3 ) 
ensemble of clouds. It would be very useful to estimate the number of absorbing clouds 
in NGC 3516. Even a continuous absorbing medium such as a wind cannot be excluded 
based on the IUE data. More importantly, high resolution data would allow accurate 
measurements of the absorbing columns; an improvement over the lower limits which can 
be established with the IUE data. If the absorbing column could be constrained to low 
values (~ 10 20 cm -2 ), then the high density clouds are probably ruled out as a location for 
the absorber: a single cloud would exceed the absorbing column, while it would produce a 
absorption feature which is far narrower than observed. 

Further analysis of archival IUE data on NGC 3516 spanning the 1978-1988 period 
is also tinder way, again using the Gaussian extraction routines. We will investigate the 
reported P Cygni profile for the absorption line variations. Also, an attempt will be 
made to measure the absorption properties of the Mg II emission line, which has not been 
reported in the literature. The Mg II line can be used to put better contraints on the 
total column density of the absorbing gas, since it samples the neutral parts of the clouds 
which are behind the illuminated faces, and which are supposed to comprise the bulk of 
the standard broad-line clouds. 

Although the X-ray spectrum did not exhibit changes in absorbing column density 
during this brief monitoring period, the X-ray flux and spectral shape did vary in a com- 
plex manner. A paper describing the simultaneous IUE and Ginga observations is in 
preparation. 
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TABLE 1 

IUE OBSERVATIONS OF NGC 3516 


Image # Int. 

time (min.) 

Date 

my 

SWP 37274 

100 

1989 Oct 7 

12.27 

SWP 37275 

250 

1989 Oct 7 

12.34 

SWP 37293 

330 

1989 Oct 9 

12.40 

LWP 16524 

65 

1989 Oct 10 

12.43 

SWP 37308 

345 

1989 Oct 12 

12.42 

LWP 16548 

60 

1989 Oct 12 

12.39 

SWP 37327 

335 

1989 Oct 16 

12.42 

LWP 16567 

60 

1989 Oct 16 

12.39 

SWP 37398 

140 

1989 Oct 19 

— 

LWP 16573 

60 

1989 Oct 19 

— 


TABLE 2 

SHORT-TERM UV VARIABILITY 
(1989 October 7 - 16) 

A 

(A) 

</> 

(mJy) 

x 2 

9 

variability® 
(% of < / >) 

1335-1360 

1.39 

4.6 

.336 

27 

1435-1470 

1.57 

7.4 

.114 

30 

1690-1730 

2.06 

7.8 

.095 

23 

1800-1842 

2.34 

10.0 

.040 

21 

1940-1970 

2.61 

15.5 

.004 

22 

o 6 

1.75 6 

1.11 

.891 

26 


a upper limit on the variability relative to < / > 
b power law index of the SWP continuum fit (F v oc u~ a ). 
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TABLE 3 

SHORT-TERM UV VARIABILITY 
(1989 October 9 - 16) 


A 

(A) 

</> 

(mJy) 

x 2 

7 

variability a 
(% of < / >) 

1335-1360 

1.32 

1.73 

.420 

20 

1435-1470 

1.49 

0.72 

.696 

13 

1690-1730 

1.97 

0.64 

.727 

10 

1800-1842 

2.27 

1.02 

.602 

8 

1940-1970 

2.51 

0.09 

.954 

6 

2655-2720 

4.68 

0.86 

.649 

10 

2960-3035 

6.07 

0.02 

.991 

7 

a b 

1.91 6 

0.45 

.800 

7 


a Upper limit on the variability relative to < / >. 

b Power law index of the IUE (SWP+LWP) continuum fit ( F v oc v ~ a ). 


TABLE 4 

EMISSION-LINE MEASUREMENTS 


ID a 

Intensity 15 

EW C 

Si IV A1397, 0 IV A1407 

16.4T1.8 

7.2T0.8 

N IV] A1484 

4.7±0.9 

1.9T0.4 

C IV A1549 

185.6T3.3 

83.7T2.4 

He II A1640 

17.8T1.0 

7.8T0.5 

0 III] A 1663 

3.2T0.5 

1.5T0.3 

N III] A1750 

3.6T0.5 

1.7T0.3 

C III] A1909, Si III] A 1892 

35.6T0.7 

16.9T0.4 

Mg II A2798 

69.0T2.2 

35.2T1.3 

0 III A3133 

<3.5±1.7 

<1.7±0.9 


a Identification of emission fine and rest wavelength. 
b Observed flux in 10 -14 ergs cm -2 s -1 . 
c Observed equivalent width in A. 
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TABLE 5 

GALACTIC ABSORPTION-LINE MEASUREMENTS 


ID 

^min (A) 

EW (A) 

0 I A1302, Si II A1304 

1302±2 

0.7T0.2 

C II A1335 

1334±1 

0.7±0.2 

Fe II A2586 

2587±3 

1.3T0.4 

Fe II A2599 

2602±3 

1.5T0.3 

Mg II AA2796, 2803 

2797±3 

1.9T0.3 


TABLE 6 

1989 ABSORPTION MEASUREMENTS 


line parameter 

Oct 7 e 

Oct 9 

Oct 12 

Oct 16 

1450 A continuum 1 

2.67±0.24 

2.18±0.20 

2.25T0.21 

2.05T0.13 

2700 A continuum a 

- 

1.91T0.12 

1.99T0.10 

1.84T0.12 

C IV blue peak flux a 

6.93 

5.76 

6.66 

6.55 

C IV red peak flux a 

6.91 

6.76 

6.74 

6.85 

C IV blue peak wavelength 15 

1542 

1542 

1541 

1546 

C IV red peak wavelength 15 

1565 

1567 

1567 

1566 

C IV minimum absorption wavelength 15 

1557 

1560 

1557 

1556 

C IV EW C 

8.2T0.4 

7.1T0.2 

8.9±0.4 

8.1±0.4 

C IV EW d 

10.4T0.3 

10.1±0.3 

10.8T0.3 

9.6T0.3 

N V A1240 EW 

7.3±0.2 

7.2T0.2 

6.8T0.2 

7.1±0.2 

Si IV A1393 EW 

2.9T0.4 

3.2T0.4 

3.6±0.4 

4.4±0.2 

Si IV A 1402 EW 

1.0T0.3 

1.8T0.3 

2.1±0.3 

1.8T0.3 


a Observed flux in 10 -14 ergs cm -2 s -1 A -1 . 
b Observed wavelength with an error of ~1 A. 

c Observed equivalent width (A) under straight line between the peaks. 
d Observed equivalent width (A) under a Gaussian fit to the narrow line component 
(FWHM~30 A). 

e SWP37275 only as SWP37274 has a cosmic ray hit at the red peak of C IV; no LWP 
spectrum obtained on this date. 
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FIGURE CAPTIONS 


FIG. 1. - Coadded spectrum of NGC 3516 obtained with IUE in 1989. Data points 
flagged with X’s are affected by cosmic ray hits and the presence of reseaux. The emission 
lines are indicated at the laboratory wavelengths redshifted with 2 = 0.009. The strong 
emission line around 1216 A is geocoronal Lya. The tick marks under the graph indicate 
probable Galactic absorption lines. The absorption features shortward of the N V and C IV 
lines are associated with the Seyfert itself. Considering the wavelengths of the absorption 
around 1400 A, these features are also thought to be associated with NGC 3516 and not 
Galactic. 

FIG. 2.- Individual spectra of NGC 3516 obtained in 1989. Note the lack of significant 
changes, except that the overall flux level on October 7 is higher than the others. 

FIG. 3.- Short-term variability of the ultraviolet continuum emission for the spectra 
obtained in 1989 October. The dashed line indicates the error-weighted average flux. The 
first set includes all the observations (see also Table 2), while the second set excludes 
October 7 (Table 3). Note the absence of significant variations at all wavelengths and the 
constancy of the power law index. 
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ABSTRACT 

We report on the first UV observations of the bright QSO HI 82 1 +643. With V = 14.2 mag and z = 0.297, 
HI 821 +643 is the second brightest object in the sky at z > 0.1. The WE data are combined with new optical 
spectroscopy, and existing infrared and X-ray data, to reveal a strong optical/UV “big bump,” which con- 
tinues past the Lyman limit in the rest frame of the QSO. A possible turnover at the high-frequency side of 
the UV continuum constrains fits of a thin accretion disk model to a large black hole mass (M =* 3 x 10 9 M G ) 
and high accretion rate ( M ~ 19 M Q yr _1 ), but a small disk size {R ou JR in — 12). The shape of the UV contin- 
uum was found to be variable, with a hardening of the spectrum when the source was brighter. Because of its 
location only 3° from the ecliptic pole, HI 821 +643 will be an important object for simultaneous UV and soft 
X-ray monitoring to test for a common origin of the UV bump and soft X-ray excess. 

Subject headings: quasars — spectrophotometry — ultraviolet: spectra 


1. INTRODUCTION 

Twenty years after the discovery of QSOs, one of the bright- 
est objects of this class had yet to be noticed. Known as 
HI 821 +643, it was identified by Pravdo & Marshall (1984) as 
the counterpart of an X-ray source first detected by HE AO 1 
and then by Einstein . With V = 14.2 mag, it is the brightest 
object after 3C 273 at z > 0.1. The brighter members of a class 
often contribute disproportionately to our knowledge of their 
physics, because they can be studied well at many wavelengths. 
The intermediate redshift of HI 821 +643 (z = 0.297) makes it 
well-suited to a study of the intrinsic ultraviolet continuum, as 
higher redshift objects suffer from absorption by intergalactic 
H I clouds. Moreover, its location only 3° from the ecliptic pole 
[the position is a(1950) = 18 h 21 m 41!67, <5(1950) = 64°19'0:8] 
makes HI 821 +643 an important object for monitoring prog- 
rams done with satellites. In particular, it is an ideal object to 
test for a common origin of the UV big bump and the soft 
X-ray excess, an association which has yet to be proven. 

In this paper, we report new ultraviolet and optical spectro- 
photometry which we combine with existing infrared and 
X-ray data to study the continuum energy distribution of this 
QSO. The combined optical and WE spectra show a strong 
UV bump. The moderate redshift allows the detection of flux 
past the Lyman limit in the rest frame, while allowing for 
optimal coverage of the UV bump within the WE bandpass. 
We present fits to a thin disk model, which require large mass 
and accretion rate for this high luminosity QSO, but also indi- 
cate a small outer radius for the disk. 

The strong soft X-ray excess, which was attributed to 
H1821+643 in the HEAO A-2 data of 1977, is absent in all 

1 Guest Observer with the International Ultraviolet Explorer Satellite. 

2 Staff member, IUE Observatory, NASA/GSFC. Postal address: Code 
668.1, NASA/GSFC, Greenbelt, MD 20771. 

3 Presidential Young Investigator. 


subsequent X-ray observations, which resemble the 
“ canonical ” hard power law. Although we take this variable 
soft excess at face value for the purpose of this paper, we also 
discuss other soft X-ray sources in the vicinity of the QSO 
which could have confused the HEAO A-2 low-energy detec- 
tor. 

2. OBSERVATIONS 
2.1. The X-Ray Identification 

When H1821+643 was identified by the Einstein imaging 
proportional counter (IPC), comparison with HEAO A-2 data 
obtained 3 yr earlier indicated that a factor of 4 decrease in the 
soft X-ray flux had occurred between 1977 and 1980 (Pravdo & 
Marshall 1984). The change could be described as the disap- 
pearance of the soft X-ray excess, which is clearly present in the 
HEAO A-2 low energy detector (LED), and not evident in the 
Einstein IPC. Subsequent observations by EX OS AT in 1984 
and 1985 (Warwick, Barstow, & Yaqoob 1989) also did not 
detect a soft X-ray excess, leading these authors to consider the 
possibility that the HEAO A-2 LED detector had been con- 
fused by another soft X-ray source in the field, possibly the 
central star of the planetary nebula Kl-16 (Grauer & Bond 
1984) which is 100" from the QSO. But having estimated the 
flux from Kl-16 in the 0.2-2 keV bandpass of the LED, 
Warwick et al. (1989) concluded that it could not have contrib- 
uted significantly to the observed soft X-ray excess. 

We can further address this issue with the Einstein IPC data 
since the IPC bandpass includes that of the HEAO A-2 LED. 
Moreover, the IPC data were taken as part of a survey of the 
north ecliptic pole region, which can be searched for other soft 
X-ray sources in the vicinity of the LED error box. There were 
four overlapping fields, all of which included HI 82 1+643, 
although in three of these the QSO fell just outside the support 
ribs of the detector, hampering detailed quantitative analysis. 
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The X-ray position and pointlike nature of the source in the 
IPC map are consistent with identification with the QSO and 
not the planetary nebula, although it cannot be excluded that 
Kl-16 might make a small contribution to the IPC X-ray 
source. More worrisome for the interpretation of the soft 
excess in the LED is another X-ray source 36' to the northeast 
of the QSO, which is of comparable intensity. We identified 
this source with the 8.0 mag K0 star, SAO 17878. The error 
boxes of the several HEAO-J detectors, in addition to the 
optical positions of possible confusing sources, are shown in 
Figure 1. The X-ray flux of SAO 17878 in the 0. 1-3.5 keV band 
as observed with the IPC is 9 x 10" 12 ergs cm -2 s~ l , about 
half that of the QSO at the time of the Einstein observations. In 
the absence of detailed knowledge about its variability, it 
cannot be ruled out that this star contaminated the LED error 
box (Nugent et al. 1983), on whose edge it lies. However, the 
X-ray flux from this star would have to have been at least 4 
times larger than the Einstein value for it to have made a 
substantial contribution to the LED flux. The only other X-ray 
source close to the LED error box is a 5.0 mag F5 star, SAO 
17828, with an IPC flux of about 1 x 10" 12 ergs cm" 2 s” 1 . 
This source is an even less likely contributor. In summary, we 
can be certain that the hard X-rays arose in the QSO, as indi- 
cated by the IPC source and the small HE AO A-I error box of 
0.05 deg 2 (Wood et al. 1984). We can also be certain that the 
soft excess was not present during the Einstein or EXOSAT 
observations. Whether the soft excess observed by the HEAO 
A-2 LED arose exclusively in the QSO cannot be determined, 
although this appears highly likely. 



RA 


Fig. 1. — Error boxes of the HEAO observations around H1821+643 
(indicated by the filled circle labeled QSO). The HEAO A-2 MED-HED error 
box (Marshall et al. 1979) is indicated by the solid lines; the HEAO A-2 LED 
error box (Nugent et al. 1983) is indicated by the dashed lines; the HEAO A-l 
error box (Wood et al. 1984) is indicated by the dot-dashed lines. The star 
indicates the position of SAO 17878, which could contribute some of the 
detected flux in the HEAO A-2 LED band. The plus sign indicates the position 
of the weak IPC source SAO 17828. The IPC error circles for the QSO and 
SAO stars have a typical diameter of 1' and are too small to be drawn in this 
figure. All positions are in 1950 coordinates. 


TABLE 1 

Einstein IPC Spectrum of HI 821 -1-643 


Date 

1980 May 25 

Exposure time (s) 

960.0 

IPC counts s-‘ (0. 1-3.9 keV) 

0.508 ±0.023 

F x (0. 1-3.5 keV)* 

1.74 x 10' 11 

L x (0. 1-3.5 keV) b 

9.6 x 10 45 

A c 

4.23 x 10" J 

a (68% confidence) 

0.60 ±0.34 

(90% confidence) 

(+0.62, -0.45) 

N h (cm " 2 ) (68% confidence) 

2.0 (+1.8, -1.5) x 10 JO 

(90% confidence) 

<5.0 x 10 JO 

(Per d o f ) 

1.34 


* Flux (ergs cm' 2 s~ ! ) in observed energy band. 
b Luminosity (ergs s~ *) in intrinsic energy band, corrected 
for absorption and assuming a = 0.6, H 0 = 50 km s _1 
Mpc~ \and q 0 — 0. 

c Normalization constant. 

2.2. The X-Ray Spectrum 

We performed a standard spectral analysis on the single 
Einstein IPC observation in which the QSO fell inside the 
support ribs. The results are listed in Table 1. Power-law fits 
yielded a best-fit energy index a = 0.6, with 68% confidence 
limits of 0.26-0.94, and 90% confidence limits of 0.15-1.12. The 
90% upper limit on the column density N H is 5 x 10 20 cm" 2 , 
which is consistent with the 21 cm value of 4.1 x 10 20 cm -2 
(Stark et al. 1990). Hard X-ray spectral measurements by 
EXOSAT (Warwick et al. 1989) and Ginga (Turner et al. 1989; 
Kii et al. 1991) are consistent with an extrapolation of the IPC 
data, both in spectral index and normalization. Therefore, it 
seems likely that the soft excess observed by the HEAO A-2 
LED is a separate component, which varies independently of 
the hard power law. 

2.3. The Ultraviolet Observations 

In 1987 we obtained the first UV spectra of H1821+643 
with the WE satellite. The strong and broad Lyman-a emis- 
sion line was saturated in the 300 minute exposure of the SWP. 
A subsequent shorter exposure (120 minutes) resulted in 
optimal detection in the Lyman-a line. All the WE observa- 
tions are listed in Table 2. The images were obtained with the 
large aperture in the low-dispersion mode, resulting in a spec- 
tral resolution of about 5 A and about 8 A for the SWP and 
LWP, respectively. 

The spectra were obtained from the line-by-line images with 
the Gaussian extraction routines at the Goddard Regional 
Data Analysis Facility to minimize the noise while maintaining 
accurate flux levels. We corrected the spectra for the effects of 
interstellar reddening, using the extinction curve of Seaton 
(1979) and E(B- V) = 0.085 mag. The color excess is based on 


TABLE 2 

WE Observations of H1821 +643 


Image 

Int. time 
(minutes) 

Date 

SWP 31431 

300 

1987 Jul 29 

LWP 11294 

130 

1987 Jul 29 

SWP 3 1523 

120 

1987 Aug 10 

SWP 35516 

275 

1989 Feb 9 

LWP 14996 

95 

1989 Feb 9 

SWP 36396 

300 

1989 Jun 4 

LWP 15659 

120 

1989 Jun 4 
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Fig. 2. — (a) Co-added SWP spectrum of H182I +643 obtained with IUE , corrected for the estimated interstellar extinction of E{B-V) = 0.085 mag. The 
Lyman-oc emission is patched in from the unsaturated image, SWP 31523. (6) LWP spectrum. The horizontal lines denote the wavelengths selected as the continuum 
and used in Figs. 4 and 6. The vertical tick marks indicate the detected interstellar absorption lines. 


the neutral hydrogen column density in the direction of 
H1821 +643 (N„ = 4.1 x 10 2 ° cm' 2 . Stark et al. 1990) and the 
Nn/EiB— V) ratio of 4.8 x 10 21 cm -2 mag -1 (Bohlin, Savage, 
& Drake 1978). We did not correct for any internal reddening, 
as there is now growing evidence that the intrinsic extinction in 
QSOs is very small (Sun & Malkan 1989a, and references 
therein). We deleted the data points which were flagged 
because of the presence of a reseau. The co-added spectrum is 
shown in Figure 2. The overall continuum level decreases 
slightly with increasing frequency. The flux level is sensitive to 
the correction for interstellar extinction, but even if the spec- 
trum is overcorrected for reddening [up to £(£ — K) = 0.2 
mag], a downward trend with shorter wavelengths remains. 
Despite the decline in sensitivity below 1200 A, there is strong 
emission short ward of the Lyman limit, which is redshifted to 
11 82 A (Fig. 3). 

As listed in Table 2, there are three series of observations 
with IUE y each resulting in at least one SWP and LWP image: 
the first in 1987, the second in 1989 February, and the third in 
1989 June. Two SWP spectra were obtained in 1987: a long 
exposure with optimal detection of the continuum level and a 
short exposure with optimal detection of the Lyman-a emis- 
sion line. In order to study the continuum, the spectra were 
binned into about 10 data points per camera, avoiding the 


major emission lines and the Galactic absorption lines. If we 
compare these three series of observations, no significant 
changes could be found between the first and second series, 
while the spectral shape did change significantly between the 
second and third series. As shown in Figure 4, the emission 
between 2000 and 3000 A decreased in intensity, while the 
short wavelength side increased by about 10%, resulting in a 
hardening of the spectrum between 1989 February and June. 

2.4. Optical Spectroscopy 

Optical CCD spectra of HI 82 1+643 covering the wave- 
length range 3138-10046 A were obtained on 1988 October 3 
UT using four different grating settings of the UV Schmidt 
spectrograph (Miller & Stone 1987) at the Cassegrain focus of 
the Shane 3 m reflector at Lick Observatory. The spectra were 
obtained through a long slit of width 8", oriented along the 
parallactic angle (Filippenko 1982). The airmass was 1.3, con- 
ditions were photometric, and the seeing was about 1 \ Typical 
integration times per setting were 300-500 s. The spectral 
resolution was 5-8 A below about 7200 A, and 10-16 A in the 
near-IR region. 

Standard techniques were used to extract one-dimensional 
spectra from the CCD data. These included subtraction of the 
bias level, division by a flat field, removal of cosmic rays, cor- 
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X (A) 

Fig. 3. — Expanded view of Fig. 2a shows the presence of substantial flux 
below the Lyman limit at 1182 A (indicated by the vertical tick-mark). We 
estimate that the covering fraction at the Lyman limit must be less than 0.09. 
The strong emission line around 1216 A is geocoronal Lyman-oc. The data 
points around 1190 A are deleted because of the presence of reseau. 


rection for geometric distortion and misalignments, and sub- 
traction of the background sky. Several stars from the lists of 
Oke & Gunn (1983) and Massey et al. (1988) were used to flux 
calibrate the spectra, as well as to remove telluric absorption 
bands at wavelengths longer than 6200 A. Excellent agreement 
in absolute fluxes and spectral shapes was found in overlap- 
ping spectral regions. However, the wavelength scale is uncer- 
tain by about 4 A, since the observations were obtained 
through a wide slit. The redshift derived from the [O in] lines is 
0.2972+0.0006. 

The overall optical spectrum of HI 821 +643, corrected for 
extinction [£(£ — V) — 0.085 mag], is shown in Figure 5. The 
near-UV flux joins that of the WE spectrum very well, as can 
been seen by comparison of Figures 2 and 5. 



Fig. 4.— Change in the UV continuum shape between 1989 February 9 
(indicated by filled circles) and 1989 June 4 (indicated by plus signs). The flux 
increased by about 10% at the shortest wavelengths while the flux decreased at 
the long-wavelength side, corresponding to a hardening of the spectrum. 


2.5. Optical and Infrared Photometry 

The UBV photometry was obtained by S. Vrtilek with the 
0.9 m telescope at Kitt Peak National Observatory on 5 nights 
centered around 1988 October 10. Variability in the V band 
was less than 0.03 mag from night to night, and less than 0.04 
mag and 0.2 mag for B — V and U — B , respectively. The 
average magnitudes are listed in Table 3. Near infrared magni- 
tudes from the Infrared Telescope Facility (IRTF) in the J (1.25 
pm\ H (1.65 //m), K (2.2 //m), and L' (3.8 /mi) bands were 
provided by M. Elvis, and are listed in Table 3. In the far 
infrared, IRAS detections in all four bands were published by 
Neugebauer et al. (1986), and are also listed in Table 3. 

The overall continuum distribution derived from all the 
above data and shifted to the rest-frame is shown in Figure 6. 
Throughout this article X and v denote observed wavelength 
and frequency, while X 0 and v 0 denote the wavelength and 
frequency in the rest-frame of the QSO. 

2.6. Radio Observations 

VLA observations at 6 cm with a partial array of six dishes 
and an exposure time of 20 minutes yielded an upper limit of 
1 5 mJy, establishing this QSO as radio quiet. 

3. RESULTS 

3. 1 . Emission Shortward of the Lyman limit 

The redshift of HI 821 +643 (0.297) is high enough to shift 
the Lyman limit to 1182 A, into the wavelength range of IUE, 
but low enough to avoid absorption by intervening H i clouds. 
Substantial flux is detected below 1182 A in all our SWP 
spectra (Fig. 2). An upper limit to the covered fraction of the 
continuum source as derived from the flux levels below and 
above 1 182 A is estimated at 0.09; a small value consistent with 
the absence of X-ray absorption in this and other bright 
AGNs. This small covering fraction of less than 0.1 is also 
expected in a high luminosity QSO such as HI 821 +643 on the 
basis of photon counting arguments concerning the photoion- 
ization of clouds on the BLR. Rest-frame Lyman continuum 
absorption has been reported (Kinney et al. 1985) in less than 
10% of all AGNs. Another bright QSO, 3C273, is a possible 


TABLE 3 

Photometry of HI 821 +643 


Band Magnitude Flux (mJy) 


U* 13.30 9.04 

B* 14.30 8.12 

V 1 14.26 7.55 

J b 13.09 8.8 

H b 12.19 13.0 

K b 11.12 22.1 

L' b 8.95 60.5 

12//m c 236 

25/rni e 373 

60/rni c 953 

\00fim e 2164 

6cm <15 


B Magnitude averaged over 5 nights, obtained with 
the automated filter photometer on the 0.9 m telescope at 
Kitt Peak National Observatory around 1988 October 
10 (S. Vrtilek 1988, private communication). 

b IRTF instrumental magnitude, obtained with 6" dia- 
meter aperture and the TnSb detector (RC2) on 1988 
April 25 (M. Elvis 1988, private communication). 

* IRAS fluxes (Neugebauer etal. 1986). 
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Fig. 5. — Optical spectrum of H 1821 + 643 obtained at Lick Observatory on 1988 October 3 under photometric conditions. The interstellar extinction correction 
for E(B — V) = 0.085 mag has been applied. The continuum appears bumpy due to the contribution of blended Fe n lines (e.g., around 3800, 4200, 5900, and 6700 A). 


cn 



Fig. 6. — Overall continuum spectrum of HI 821 -1-643 plotted in the QSO 
rest frame, with the co-added WE data from 1987. Note the turnover of the 
UV bump and the large amplitude variation between the X-ray observations. 
The power-law plus thin accretion disk (solid line) fit with M — 3 x 10 9 M 0 , 
tii = \9 M q yr _! , and R^JR^ — 12, deviates r between log v 0 = 14.8 and 
log v o — 1 5.2 due to the Balmer continuum and Fe u line emission. 


exception, as a Lyman discontinuity has been reported in com- 
bined I U Ef Voyager data (Reichert et al. 1988). By contrast, the 
continuum spectrum of HI 82 1+643 seems to turn over 
smoothly well longward of the Lyman limit, which argues that 
the thermal peak, rather than any Lyman limit absorption, is 
responsible for the drop in flux level at the high-frequency side 
of the IUE spectrum. However, one must be aware of the 
possibility that Lyman continuum absorption is an intrinsic 
property of accretion disks in AGNs (Sun & Malkan 1987; Sun 
& Malkan 1989b; Scott & O’Dell 1989), and that a suitable 
smearing of the Lyman edge due to strong Doppler and gravi- 
tational effects in the inner disk might cause some of the turn- 
over which is observed. 

3.2. Continuum Energy Distribution 

None of the observations in different wavebands were simul- 
taneous; nevertheless, the infrared through UV data are very 
continuous. There is a steep rise in the flux level at optical 
frequencies, and the IUE data points indicate a strong UV 
bump. The prominence of the UV bump is similar to the one 
detected in the optical/UV bump prototype PG 1211 + 143, 
but it could be observed to higher rest-frame frequencies due to 
the higher redshift (z = 0.297 vs. z = 0.085 for PG 1211 + 143). 
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In all the spectra except the one obtained in 1989 June, there is 
evidence for a turnover at the high-frequency side (around 
log v 0 = 15.3). 

Because of the 5 pm feature displayed by the IR data points, 
the power law, if any, is not well defined in the infrared. By 
forcing a fit to the 60 //m and 1 pm data points, we find 
a = 1.16 (where F v oc v~ a ). Because there is a strong correlation 
between the IR and hard X-rays in QSOs (Edelson Sc Malkan 
1986; Kriss 1988), a single power law might be extrapolated 
from log v 0 % 12 to log v 0 » 18. The turn-up from the overall 
power law of the hard X-ray data (both for the Einstein data 
and the 2-10 keV HEAO 1 data) is typical for the high- 
frequency side of the continuum energy distribution in QSOs 
(Elvis 1988). Considering the large amplitude of variability in 
the X-rays, it is impossible to estimate the soft X-ray flux level 
at the time of the IUE observations, which would be important 
for models of the soft X-ray excess. 

3.3. Emission Lines 

Information on continuum emission in the unobservable 
EUV gap can be derived from the UV and optical emission 
lines, which are powered by photons with EUV energies. To 
this end, we measured the major emission lines in the UV and 
optical spectra. No significant differences in line fluxes were 
found from one IUE observation to the other; therefore, all the 
IUE spectra were coadded and dereddened. The measured line 
properties are listed in Table 4 and the lines are identified in 
Figures 2 and 5. The Lyman-a emission line could only be 


TABLE 4 


Emission-Line Measurements of HI 821 + 643 


ID 1 

Intensity b 

Flux c 

EW' 

(A) 

Cm 2977 

<35,0 

<20.3 

<5.0 

O vi 21034* 

137.1 

83.4 

16.8 

Lya 21216, (Nv 21240) 

963.5 

626.9 

156.0 

Si iv 21397,0 iv 21407 

49.8 

32.8 

10.0 

Civ 21549 

504.8 

305.9 

111.4 

He ii 21640, (Om] 21663) 

161.9 

91.7 

36.8 

Cm] 21909 

115.2 

78.1 

30.4 

Mg ii 22798 

222.6 

186.5 

98.7 

Om 23133 

1.8 

1.6 

1.0 

[Nev] 23345 

2.1 

1.8 

1.3 

[Ne v] 23426 

4.3 

3.8 

2.8 

[O n] 23727 

5.1 

4.6 

4.1 

[Ne m] 23869 

4.5 

4.2 

4.7 

H£ 23889, He 1 23888 

2.4 

2.2 

2.5 

[Ne hi] 23968, H* 23970 

3.6 

3.3 

3.6 

H<5 24102 

11.0 

10.3 

13.2 

Hy 24340, [O m] 24363 

29.0 

27.8 

39.3 

He ii 24686 

0.6 

0.6 

1.0 

H/? 24861 

100.0 

100.0 

172.6 

[O hi] 24959 

8.4 

8.3 

10.4 

[O HI] 25007 

28.6 

28.4 

40.3 

He 1 25876 

9.8 

10.2 

24.0 

[O i] 26300 

1.6 

1.7 

4.0 

[O i] 26363, ([Fex] 26375) 

0.9 

0.9 

2.0 

Ha 26563, ([Nil] 226548,83) 

333.3 

352.3 

895.1 

[S ii] 267 16,30 (He 1 26678) 

4.8 

5.1 

9,9 


* Identification of emission line and rest wavelength; the lines in parenth- 
eses denote only minor contributions. 

b Intensity corrected for Galactic reddening, relative to with 
/(H/0 = 1.05 x 10~ 12 ergscm -2 s -1 . 

* Observed flux relative to H/7 with F(H/7) = 8.63 x 10 _ 13 ergs cm -2 s" L 

* Observed equivalent width. 

* Narrow component only (FWOI —40 A). 


measured in SWP 31523, as the other SWP spectra were satu- 
rated at this line. The major source of error in the line measure- 
ments is the uncertainty in the continuum level, which is 
contaminated by the broad wings of the emission lines, Galac- 
tic absorption lines, and blended Fe n line and Balmer contin- 
uum emission. For the O vi 21034 emission line only the 
narrow component was measured (with a full width at zero 
intensity, FWOI, of about 40 A), as it was unclear to where the 
broad wings extend (possibly from 1280 to 1380 A, encompass- 
ing the Galactic absorption feature around 1300 A). This 
uncertainty around the O vi line also complicated the mea- 
surement of the upper limit to C nr 2977 at its redshifted wave- 
length of 1267 A, as well as the measurement of the O i 21302, 
Si ii 21304 absorption feature. 

These UV observations are the first to be conducted for 
H 1821 +643. There were, however, previous optical observa- 
tions with emission line measurements from 1981 (Pravdo Sc 
Marshall 1984). The equivalent widths listed by these authors 
for Mg ii, Hy and Hfi are smaller (by about 30%) than our 
measurements, and larger (by 15%) for [O hi] 25007. Com- 
pared to a sample of 20 intermediate redshift QSOs (z < 0.84, 
Kinney et al. 1985) the EW of Lyman-a for H1821+643 is 
typical, being close to the sample’s average. Kwan Sc Krolik 
(1981) compiled the average and range of permitted line ratios 
for a large sample of QSOs. No single line ratio observed in 
HI 821 +643 is outside this sample’s range. Compared to the 
mean, the observed ratios of Lya/H/?, C iv 21549/Lya, and 
Mg ii 22798/H/? in HI 821 +643 are large (by a factor of about 
1.5), while the C in] 21909/C iv 21549 ratio is half the mean 
value. 

Considering the overall energy distribution of HI 821 +643, 
it would be of interest to determine where the spectrum peaks: 
at IUE frequencies as suggested by the detected turnover, at 
soft X-ray energies as indicated by the HEAO I data, or in the 
unobservable EUV gap. Krolik & Kallman (1988) investigated 
this question theoretically by modeling the effects of the several 
different accretion disk spectra [i.e., bare power law, power law 
plus disk spectrum peaking at 10 eV (log v 0 ~ 15.4), and one 
peaking at 80 eV (log v 0 — 16.3)] on the emission line forma- 
tion. Although the ionization parameter E is far more impor- 
tant than the shape of the ionizing spectrum in determining the 
line emission, there are line ratios which are relatively insensi- 
tive to the ionization parameter. 

The observed line ratios in HI 821 +643 are consistent with 
a power law plus disk spectrum peaking at 10 eV, i.e., in the 
UV, within the context of the Krolik Sc Kallman (1988) model. 
All ratios are within the model range (with the pressure as a 
free parameter between 6.3 x 10 13 and 20 x 10 13 K cm' 3 ) for 
the high ionization case (S = 2.0), except for C m] 21909/C iv 
21549 (the one ratio that is substantially smaller than the 
average value for Kwan Sc Krolik’s 1981 sample of AGNs). A 
high E might be expected for an ionizing spectrum peaking at 
low energies in order to produce enough high-ionization 
species. 

3.4. Absorption Lines 

The location of HI 821 +643 at (/, b) = (94°, 27°) and the 
moderate column density of 4.1 x 10 20 cm” 2 are ideal for 
study of interstellar absorption. The Galactic longitude of 
HI 82 1+643 permits optimal velocity separation of distant 
matter from local matter, while the low-latitude permits the 
viewing of distant matter at large height above the Galactic 
plane. Therefore, future high-resolution studies of the absorp- 
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TABLE 5 

Absorption-Line Measurements of H 1 82 1 + 643 


EW 

id (A) 


S ii A1 25 1 , 1254, 1260, Si ii A 1 260 0.8 

O l ^ 1 302, Si ll 21304 2.6 

Civ >11548,51 <0.7 

Ain >11671 0.6 

Fe ii >12343 2.8 

Fc ii 22374,2382 3.5 

Fen 22586 1.6 

Fen 22599 1.8 

Mg ii 22796,2803 3.9 


tion features in HI 821 +643 could investigate the properties of 
the interstellar medium at large galactocentric distances. We 
measured the absorption lines listed in Table 5 and indicated 
by tick marks in Figure 2. The lines are from low-ionization 
species, except for C iv for which only an upper limit is 
obtained due to its location on the blue wing of Lyman-a. The 
C iv A 1548, 1551 absorption is thought to arise in the Galactic 
halo (Savage 1988). We did not detect any absorption lines in 
the optical spectrum. An extensive compilation of EWs of 
absorption lines arising at higher Galactic latitude ( b ~ 60°, 
Savage 1988) gives lower values than ours (typically by a factor 
of about 2) for the low-ionization species except S ii. Si ii and 
AI n. Other studies of absorption features arising at lower 
latitudes [(/, b) = (96°, 12°), Pettini et al. 1982] show the same 
strong absorption in Mg n 22796,2803 and Fe n 22599 as we 
measu red in H 1 82 1 + 643. 

In H 1821 +643, the spectral region between Lyman-a and 
C iv 21549 covers the redshift range 0.04-0.27 for detection of 
C iv absorption lines arising between the QSO and our 
Galaxy. In this sensitive part of the SWP camera no such 
absorption features were found above the detection threshold 
of EW — 0.3 A. The spectral region shortward of the Lyman a 
emission line allows the detection of low-redshift Lyman a 
forest lines. Although in the co-added SWP spectrum there are 
several possible absorption features in this region (e.g., at 
2 ~ 1470 and 1490 A), none were present persistently in all four 
SWP spectra. The absence of intergalactic Lyman-a and C iv 
lines is expected for the low redshift of HI 821 +643 (Bergeron 
1988). 

4. ACCRETION DISK MODELS 

Through studies of variability at different wavelengths (Cutri 
et al. 1985), and continuum energy distributions between IR, 
optical, UV, and X-rays (Malkan & Sargent 1982), it has been 
established that the optical/UV bump is a separate feature 
from the underlying power law and is now usually attributed 
to an accretion disk. 

In the original model of Shakura & Sunyaev (1973), a sum- 
mation of Planckian emission arises from the geometrically 
thin (but optically thick) accretion disk. At the high-frequency 
side a turnover occurs, corresponding to that blackbody curve 
with the maximal temperature at the inner disk. There are two 
free parameters in this model (in addition to the inclination i, 
and the outer radius of the disk, which affects the low-energy 
side): the mass M of the central black hole and the accretion 
rate M. In addition, the optical continuum generally deviates 
from the overall power law around 3000 A. This emission is 
attributed to Balmer continuum and blended Fe ii line emis- 
sion (Wills, Netzer, & Wills 1985). 


Several modifications to the original model have been sug- 
gested: inclusion of opacity (e.g., electron scattering and 
Comptonization) and inclination effects (Czerny & Elvis 1987; 
Wandel & Petrosian 1988). The Comptonization hardens the 
emitted spectrum without increasing the energy requirements, 
and therefore keeps the luminosity sub-Eddington. One of the 
inclination effects is relativistic focussing of the innermost rays 
along the disk (Cunningham 1975), resulting in harder spectra 
for disks which are viewed edge-on. The effects of scattering 
around the disk have been modeled using stellar atmospheres 
(Sun & Malkan 1987), resulting in a strong Lyman limit 
absorption edge for disks observed face-on. However, these 
authors indicate that blackbody disk models fit the observa- 
tions better than the models using stellar atmospheres (Sun & 
Malkan 1989a). 

Although alternatives to the thin accretion disk models are 
available (e.g., the thick torus model by Madau 1988), thin 
disks are widely used as these have been shown to fit observa- 
tions of QSOs and Seyfert 1 galaxies successfully (Malkan & 
Sargent 1982). We calculated disk models for a nonrotating 
black hole (a/M = 0, the Newtonian case), with the standard 
assumption that the emerging radiation can be described 
locally by a blackbody spectrum. The temperature T at radius 
R is given by 


T(R) = 



(1) 


where R in is the innermost stable orbit ( R in = 3R Sch , Shakura 
& Sunyaev 1973). We did not include relativistic effects, since 
the disk around a nonrotating black hole does not reach close 
to the black hole (unlike the case of the rotating black hole, 
where the disk reaches much closer) and the relativistic contri- 
butions are, therefore, relatively small. The emergent flux is 
calculated after fixing the inclination angle such that 
cos i = 0.5 and assuming a distance D given by 



for q 0 = 0 and H 0 = 50 km s“ 1 Mpc~ *. We did not include a 
fitting component to account for the Balmer continuum and 
blended Fe ii line emission, and, therefore, our fits deviate from 
the observed data points between log v 0 ~ 14.8 and log v 0 ~ 
15.2. 

The best-fit (by eye) to the 1987 IV E data, as shown in 
Figure 6, with a thin accretion disk model added to the power 
law, is satisfactory to explain the optical/UV spectrum with the 
following fitting parameters: M = (3 + 0.5) x 10 9 M Qy M = 
19 ± 3 M Q yr- 1 , R„JR in = 12( + 5, -3), (L/L Md = 0.16). The 
fit is constrained at both the high-energy side because of the 
turnover of the high-frequency IUE data points and at the 
low-energy side due to the steep rise of optical data points. The 
major sources of uncertainty in the fitting parameters are (1) 
the nonsimultaneity of the optical and UV data points, (2) the 
contribution of the Balmer continuum and the Fe 11 lines, and 
(3) the correction for interstellar reddening. In this model, the 
mass of the central object determines mainly the frequency of 
the peak, while the accretion rate determines mostly the nor- 
malization. Once these two parameters have been determined 
by the turnover in the UV and the strength of the optical/UV 
bump, respectively, the width of the bump is determined by the 
extent of the disk, i.e., the ratio R ou JR in . The lower limit to 
Roui/^in rules out the possibility of a single blackbody fit, which 
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would result in a far too narrow bump. If the reality of the UV 
turnover seen around log v 0 = 15.3 is doubted, then the bump 
could peak at higher frequencies, outside the range of IUE. In 
that case, M would be smaller and the bump would be wider 
requiring a larger disk (possibly increasing R ou JR in to about 
30). Such a result is within the range of acceptable fits to the 
1989 June IUE data, as the short-wavelength turnover is not 
apparent then. To estimate the uncertainty due to the inter- 
stellar reddening correction, we found that E(B— V) = 0.1 mag 
(an increase by about 20%), would increase M by about 10%, 
while M and the size of the disk would remain unchanged. The 
use of a color excess of E(B— V) — 0.12 mag or higher would 
result in a dereddened spectrum without a turnover, again 
increasing the uncertainty in M and R oux /R in . However, such a 
large color excess would correspond to a hydrogen column 
density of 5.8 x 10 2 ° cm -2 , far above the values measured near 
HI 821 + 643. 

The novel result of this analysis is the relatively small value 
of R out /R {n (9-17). Bechtold et al. (1987) applied a similar fitting 
model, to the spectrum of PG 1211 + 143 and found a best fit 
Koui/^in of at least 200, considerably larger than we found in 
HI 82 1 +643. It is of interest to note that the unexpectedly 
small outer radius in HI 821 +643 could be identified with the 
maximum or critical radius (R cri( ) in the so-called two- 
temperature accretion flow model of Begelman, Sikora, & Rees 
(1987). These authors studied the ion torus around a central 
black hole, and in particular, the cooling effects of electron- 
positron pair production on such a torus. Due to these cooling 
effects the torus collapses to a geometrically thin, optically 
thick annulus in the innermost region, defined by R < R cri( ~ 
30 R The thermal UV and possibly the soft X-rays would 
come from this optically thick region, which would in fact have 
an outer radius remarkably close to our fit value R out = 36 R Sch 
(R ou JR in = 12). In this model, due to the bistability of the 
annulus, smooth changes in M would cause abrupt changes in 
the flow structure, which in turn would be responsible for spec- 
tral variations. In a stable disk, changes in M would affect the 
emergent spectrum only on the viscosity time scale t vh 
(Bechtold et al. 1987); for an a disk with the parameters 
inferred for HI 821 +643 above, f vis is greater than 6000 yr for 
a < 1.0 (Frank, King, & Raine 1985). 

A successful fit of both the optical/UV bump and the HE AO 
I soft X-ray excess is practically impossible to achieve, con- 
sidering the decline of the IUE data points, and the high and 
fairly constant flux level for the 0.5-2 keV HEAO 1 observa- 
tions. If such a fit is forced (ignoring the nonsimultaneity of the 
UV and X-ray data), highly super-Eddington luminosities 
(L ~ 40 L cdd ) would be produced for the Newtonian accretion 
disk models, violating the thin disk approximation for the 
inner disk (Bechtold et al. 1987), unless Compton scattering of 
the UV photons to EUV and soft X-ray energies is incorpo- 
rated (Czerny & Elvis 1987). While following the optical data 
points rather closely, such a fit exceeds the observed UV flux 
by a factor of about 2-3. Moreover, as there is no evidence that 
the UV emission and the soft X-ray excess are linked in AGNs 
(Elvis, Wilkes, & McDowell 1989), there is no convincing justi- 
fication to fit both features in one model fit for HI 821 +643 at 
this moment. This QSO would, however, be an excellent candi- 
date to study whether the UV and soft X-ray radiation have 
the same physical origin, considering the variability and 
brightness of both the UV and X-ray emission. In Figure 4 the 
change in the UV continuum over a period of 4 months is 
shown. This change in shape, where the spectral hardness is 


correlated with luminosity, is often found in AGNs which are 
well observed with IUE (Edelson, Krolik, & Pike 1990) and is 
consistent with an expected increase of temperature with accre- 
tion rate. The reported minimum doubling time of about 
10 days at X-ray energies (Snyder & Wood 1984) could be 
identified with the dynamical time scale (f dyn ) of the inner 
disk. We define f dyn as (I + z) R ma > max , where R mBX is the 
radius at which the temperature of the disk reaches a 
maximum (R ~ 8 GM/c 2 ) and = ( GM/R m J 1/2 . Thus, 

i dy „ = GM(8/c 2 ) 3/2 (1 + z) as seen by the observer, or t d ~ 1.7 
(M/10 9 M G )days. 

Sun & Malkan (1989a) presented fits of a thin accretion disk 
model to the spectra of 24 AGNs, not including HI 821 +643. 
As the inferred disk parameters are very similar in the Sch- 
warzschild metric (used by Sun & Malkan 1989a) and the 
Newtonian metric (used here) for cos i — 0.5, we can compare 
our values obtained for mass and accretion rate with the large 
sample. With Af = 3 x 10 9 Af 0 , H1821 +643 would be one of 
the most massive AGNs in the sample and the most massive of 
the AGNs with intermediate redshift (z = 0.2- 1.0). The inferred 
accretion rate of M = 19 Af 0 yr “ 1 would be one of the highest 
for an intermediate redshift QSOs. This was to be expected, 
considering the large luminosity of HI 821 +643 and the pro- 
nounced optical/UV bump in its spectrum. 

There is, however, one significant difference between our 
model and that of Sun & Malkan (1989a), and that is the use of 
the outer radius of the accretion disk. Sun & Malkan (1989a) 
do not specify their outer radius, but it can be assumed to be a 
large value. In our fitting model the outer radius is a free 
parameter, resulting typically in low inferred values for suc- 
cessful fits (R ou JR in ~ 12). A strong UV/optical bump with a 
small width, as is the case for H 1 82 1 + 643, can only be fitted in 
a satisfactory manner with a small accretion disk. While the 
turnover and the height of the UV bump determine the central 
mass and accretion rate, the extent of the disk is determined by 
the width of the bump. We maintain that only a limited disk, in 
which the range of temperatures is small and therefore the 
emergent spectrum only slightly broader than the single black- 
body curve, fits a narrow optical/UV bump, as observed in 
HI 821 +643. 

5. CONCLUSIONS 

We have obtained multiwavelength observations of 
HI 821 +643 and have modeled its continuum with an accre- 
tion disk in addition to an underlying power law. The contin- 
uum follows the power law with a = 1.16 most closely in the 
far-IR and near-IR, deviating strongly at optical and UV fre- 
quencies, and possibly at soft X-ray energies. Most of the 
energy is emitted in the UV, while the variable soft X-ray 
excess could have a similar luminosity when at maximum. Fits 
with a geometrically thin, optically thick accretion disk yield a 
central black hole mass of 3 x I0 9 M 0 and an accretion rate of 
19 M 0 yr -1 , both large values consistent with the luminosity 
of HI 821 +643 and the prominence of the optical/UV bump in 
its continuum spectrum. The accretion disk fits result in a 
surprisingly small disk size (R ou JR in = 12). Efforts to constrain 
the extreme UV spectrum by applying photoionization 
models, which predict emission line ratios, were indicative of 
an ionizing spectrum peaking at UV frequencies. 

The shape of the UV continuum spectrum was found to vary 
significantly, but the time scale of the UV variability remains to 
be determined by periodic observations. An outstanding ques- 
tion for QSOs with both a strong optical/UV bump and a soft 
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X-ray excess, of which H1821+643 is a prime example, is 
whether these components have a common origin. Simulta- 
neous observations at UV and soft X-ray frequencies could 
determine whether both features are present at the same time. 
If both emissions arise in an accretion disk, then this should be 
reflected in the continuum energy distribution and the corre- 
lated variability of the UV and soft X-ray flux. To fit both the 
UV bump and the soft X-ray excess as one feature seems prob- 
lematic in the current accretion disk models (e.g., requiring 
highly inclined disks around rapidly rotating holes, as shown 
by Sun & Malkan 1989a, or the inclusion of a strong 
Comptonized component, as calculated by Czerny & Elvis 
1987). HI 821 +643, with its strong emission in all wavelength 
bands (except radio), and its variability in the UV and soft 
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